High-harmonic transient grating spectroscopy of NO2 electronic relaxation 
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We study theoretically and experimentally the electronic relaxation of NO2 molecules excited by absorption 
of one ~ 400 nm pump photon. Semi-classical simulations based on trajectory surface hopping calculations 
are performed. They predict fast oscillations of the electronic character around the intersection of the ground 
and first excited diabatic states. An experiment based on high-order harmonic transient grating spectroscopy 
reveals dynamics occuring on the same timescale. A systematic study of the detected transient is conducted to 
investigate the possible influence of the pump intensity, pump wavelength, and rotational temperature of the 
molecules. The quantitative agreement between measured and predicted dynamics shows that, in NO2, high 
harmonic transient grating spectroscopy encodes vibrational dynamics underlying the electronic relaxation. 



I. INTRODUCTION 

Resolving in time the energy relaxation in molecules 
both at the vibrational and electronic levels is theo- 
retically and experimentally challenging. There are a 
variety of experimental techniques for probing molec- 
ular dynamics, whose utility can vary from one poly- 
atomic molecule to another one 1 ' 4 . Among the most 
commonly used pump-probe signals, we can highlight 
photoelectron energy^and angular distributions (PAD) 6 , 
electron-ion coincidence in the perturbative regime-^ 
that can provide molecular frame- PAE0H1 or in the non- 
perturbative regime like Coulomb explosion^ or above 
threshold ionization^]. All these detections are based on 
collecting charged species. An alternative experimental 
technique developped in the last decade is high harmonic 
spectroscopy, based on the high order harmonic genera- 
tion process. This method has shown a relevant sensitiv- 
ity to nuclear dynamics^HXH anc j an inherent connection 
to the symmetries of molecular orbitali^^. 

The high harmonic generation (HHG) process, which 
occurs when a molecule is submitted to a strong probe 
laser field, can be described in a first approximation by a 
three step modePSl. A valence electron first tunnels out 
of the molecule under the influence of the strong field. 
This electron is then driven away by the electric field and 
can finally recombine onto its ionic core to emit extreme 
ultraviolet radiation. This process occurs every half pe- 
riod of the probe laser, so that the resulting radiation 
is constituted of odd high-order harmonics. The ampli- 
tude, polarization and phase of the XUV-VUV emission 
encode information about the electronic configuration of 
the target molecules. 

Since this pump-probe signal is purely based on 
optical detection, its contrast can be drastically re- 



duced due to the harmonic emission from unexcited 
molecules. This can be improved by using for instance 
polarization resolved spectroscopy 2 ^ or transient grating 
spectroscopy 2 ^. In the latter case, a spatial grating of 
molecular excitation is created by optical interferences 
between two identical non-colinear pump pulses. The si- 
nusoidal spatial modulation of the pump intensity across 
the molecular beam results in a partial diffraction of the 
high harmonic emission. 

This approach has been used in the present paper to 
revisit the relaxation dynamics in NO2 molecules elec- 
tronically excited in the first optically-allowed electronic 
state 2 B2 close to the first dissociation limit (~ 3 eV, 
Fig. [I]). The dipole of this transition belongs to the B 2 
irreducible representation of the C2,; point group, namely 
the pump pulses preferentially excite NO2 molecules with 
the C2„ axis perpendicular to the pump polarization. By 
reducing the point group symmetry to C s , the l 2 Ai and 
I 2 B2 states become the 1 2 A' and 2 2 A' states respectively 
that are vibronically coupled. We call these diabatic 
states 1 and 2, respectively, as labelled in Fig. I. This 
strong coupling results in a high density of states (0.754 
levels/cm - and leads to the loss of both the elec- 
tronic and vibration/rotation signatures. This has made 
of NO2 a kind of hellish grail of molecular physicists. 
Time-resolved photoionization has revealed unassigned 
regular oscillati ons w ith a 500 fs period for a pump cen- 
tered at 400 mdSEBEZ] 

On the contrary, a recent high harmonic spectroscopy 
experiment using ~ 400 nm pump pulses has shown a few 
oscillations of ~ 100 fs characteristic time that have been 
assigned to the 2 B2- 2 Ai vibronic coupling^in qualita- 
tive agreement with theoretical prediction d 29 * 30 l In the 
present paper, we reproduce this experiment and find 
subtle but significant differences in the femtosecond tran- 




FIG. 1. Potential energy surfaces of the diabatic fundamental 
state 1 2 A' and first excited state 2 2 A' of NO2 as a function of 
the bending angle and NO strecth distance Ri . These diabatic 
states are referred to as 1 and 2 in the text, respectively. The 
second NO distance R2 is fixed to 1.19 A. 

sient with respect to previous results^. We investigate 
the influence of the pump intensity and wavelength, and 
of the rotational temperature of the molecular beam. 
Simulations based on Trajectory Surface Hopping (TSH) 
reveal that the diffracted signal from the excitation grat- 
ing directly reflects the vibrational dynamics of the first 
excited diabatic state around the conical intersection. 
The remainder of the paper is organized as follows. Cal- 
culations of the molecular dynamics are presented in Sec- 
tion II. The experimental results of high harmonic tran- 
sient grating spectroscopy are given in Section III. The- 
ory and experiments are combined to deliver an interpre- 
tation to the observed dynamics in Section IV. 



II. SIMULATIONS 

A. Background 

NO2 has received much attention because it is a sta- 
ble and small molecule for which the effect of a coni- 
cal intersection can be observed and calculated. Most 
of the theoretical studies were focused on the analysis of 
the experimental spectra and were performed in a time- 
independent approach. More recently, time-dependent 
simulations of the intramolecular nonadiabatic dynamics 
were also performecP^il. These studies were based on 
a diabatic representation of the ele ctroni c Hamiltonian 
obtained from ab initio calculation ;] 31 ! 32 ! ' or from the 
experimental spectrum 3 ^. By propagating a quantum 
wavepacket or a swarm of surface hopping trajectories, 
the time dependent populations of diabatic or adiabatic 
states were calculated and regular recurrences were ob- 
served in some of these populations. For a range of total 
energy going from ~ 13000 cm -1 up to ~ 28000 cm -1 , 
the period of these recurrences is in the range 50 - 100 fs 
with a propensity to increase when total energy increases. 



A good agreement was observed be tween quantum time- 
dependent (QTD) and TSH results^ 21221 , indicating that 
the semiclassical approach is relevant for this system. 

Since experiments are carried out around 3 eV, the 
comparison with theory requires accurate PES including 
the O+NO dissociation channel. Among the three exist- 
ing models of PESpHSfD satisfying these requirements, the 
most reliable one was built by Kurkal, Fleurat-Lessard 
and Schinke 36 (KFS) from ab initio data calculated with 
multireference configuration interaction and a quadruple 
zeta basis set. This last work provides both the adia- 
batic and diabatic representations of the ground and the 
first excited PES. The bound levels of the ground elec- 
tronic state were calculated up to ~ 7050 cm -1 and a 
good qualitative agreement with the experimental spec- 
trum was obtained. 

However the KFS diabatic representation is incomplete 
because neither the diabatic electronic coupling nor the 
nonadiabatic derivative coupling have been calculated. 
Nevertheless the square of the electronic coupling V12 can 
be extracted from the available diabatic (Ei and E 2 ) and 
adiabatic (Ex and Ea) energies by 

V 1 \ = \[(E A -E x f-(E 2 -E l f] (1) 

Let Q_ = B,\ — B,'2 be the antisymmetric stretch coor- 
dinate, Ri and R2 being the N-0 distances. When Q_ 
= 0, the molecule has the Civ symmetry, the electronic 
coupling should be equal to zero and the adiabatic ener- 
gies should be equal to the diabatic energies. If Q_ 7^ 0, 
the symmetry is reduced to Cs, the electronic coupling 
can be different from zero and the following inequality 
must be fullfiled: 

\E A -E X \>\E 2 -E l \ (2) 

Because the adiabatic and diabatic PES have been sepa- 
rately fitted, the inequality ^ is not always satisfied due 
to the small errors introduced by the fitting procedures. 
Since this occurs in regions where V12 is close to zero, we 
safely reduce it to exactly zero. The sign of V12 cannot 
be known from equation ((TJ, but at the first order, V12 
is proportional to Q _I 37 I 3S T p or this reason which is a 
particular case of a more general symmetry reason 3 ^, V12 
takes the sign of Q_. 

Finally, let us remind two important energy quantities 
of the PES. The O+NO dissociation energy of the ground 
electronic state is Do = 3.00 eV while the experimental 
value is Do = 3.11 eV^. The difference of potential en- 
ergy from the ground to the excited electronic state, both 
at equilibrium configuration, is T e = 1.31 eV, while the 
experimental value is 1.21 eV^B. 

B. Trajectory surface hopping calculations 

Although the nonadiabatic dynamics of a triatomic 
molecule can be simulated nowadays with a QTD method 
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without difficulties, we have preferred the semiclassical 
approach. The obvious benefit of using semiclassical 
methods is the possibility to study large polyatomic sys- 
tems for which quantum dynamics calculations are out 
of reach. Validating the semiclassical approach for the 
small system NO2 will pave the way to the study of 
larger systems. While TSH simulations of the dynam- 
ics of NO2 have been already successfully compared to 
QTD simulation d 32 * 33 !, the present work provides the first 
comparison of the former with experimental data. The 
TSH method used here is based on the Tully's fewest 
switches algorithm^ and some implementation details 
can be found in a previous worlP^l. 

The starting point in the phase space and the start- 
ing time of the trajectories are obtained by applying the 
classical treatment proposed by Meier and EngeP^. The 
initial phase space distribution is obtained by a two-step 
procedure. First, the initial distribution for the vibra- 
tional ground state of the electronic ground state, rep- 
resented by the Wigner function p w , is filtered by the 
resonant one-photon transition 

P(q,p)=^ ) (q,p)e- (15l2 - ftw)2/2aft2 (3) 

where D12 is the potential energy difference of the tran- 
sition, uj the central frequency of the pulse, and a the 

.2 

parameter of the pulse envelope defined by e~ . Sec- 
ond, the random starting time of each trajectory is drawn 
from the distribution function e~ 2at . The quantities 10 
and a have been defined by considering a pump pulse 
centred at 400 nm and with a 30 fs FWHM gaussian 
envelope. 

Figure [2] shows two different evolutions of the popu- 
lation P 2 of the excited diabatic state: the first one is 
calculated by applying only the first step of the initial 
procedure, and the second one is calculated with the com- 
plete initial procedure. As expected, using random start- 
ing time tends to smooth out the details of the dynam- 
ics. These calculations have been done by propagating a 
batch of 1 x 10 4 trajectories. Smaller batches of 2 x 10 3 
and 5 x 10 3 result in no observable difference in the P2 
population behavior, except that some numerical noise 
increases progressively as the number of trajectories de- 
creases. The influence of the pump wavelength was also 
investigated by computing P2 for a laser pulse centred at 
395 nm and at 392 nm. No significant difference could 
be observed with respect to the result at 400 nm shown 
111 Fig [2j All three curves were almost indistinguishable 
until t = 100 fs. 

Let us now examine the intramolecular nonadiabatic 
dynamics revealed by these calculations. Figure [3] shows 
the evolution of the populations P± and P2 on the ground 
and excited diabatic PES respectively, along with the 
evolution of the average bond angle on each PES. Addi- 
tionally, Figure [4] displays some snapshots of the trajec- 
tories swarm on both diabatic surfaces at selected time 
steps. Because the details of the dynamics are more eas- 
ily observed in the case where random initial delay is 
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FIG. 2. Temporal evolution of the P2 diabatic population 
calculated with a batch of trajectories starting at the same 
time (green) or starting at gaussian random time (red), in 
the case of a 30 fs pump pulse centred at 400 nm. 




FIG. 3. Temporal evolution of the ground (blue) and first 
excited (red) diabatic state populations (a) and average bond 
angle (b) from TSH calculations, assuming a 100% excitation 
probability. 



not applied, the snapshots where obtained from a batch 
of trajectories starting at the same time on the excited 
PES. Only 2 x 10 3 trajectories are included in order to 
avoid overloading the figure. At the beginning, Pi(t = 0) 
= and P2(t = 0) = 1. Then the swarm goes downhill 
and encounters the crossing seam at t ~ 10 fs. Because 
the swarm is still quite close to a C% v geometry, the elec- 
tronic coupling is inefficient, with few trajectories jump- 
ing on the ground PES. The swarm reaches the turning 
point at small bond angles at t ~ 20 fs, spreads along Q_ 
before returning to the crossing seam. Due to this Q_ 
dispersion, the electronic coupling is effective at this sec- 
ond crossing of the seam. Therefore, a large number of 
trajectories jumps on the ground state. The main effect 
is the drastic decrease of P2 observed between t = 30 fs 
and t = 60 fs on Figure 3 (a) . The variation of P2 is now 
controlled by the bending motion of the swarm on the 
ground PES. This swarm moves to the O-N-0 linear con- 
figuration (t — 60 fs) and then goes back towards small 
angles, crosses again the seam and thus repopulates the 
excited PES: this is the first recurrence of P2 at t ~ 100 
fs. Another period of bending occurs, indicated by a sec- 
ond recurrence which is hardly visible around ~ 170 fs. 
Then both swarms are completely scattered in the whole 
phase space and no more change of the populations can 
be seen. 

In the adiabatic representation, the dynamics is some- 
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FIG. 4. Snapshots of the trajectories swarm and contour lines 
of the diabatic PES. The blue and red dots represent the coor- 
dinates of the trajectories on the ground and excited diabatic 
PES respectively. The contour lines are separated by 0.5 eV 
and the blue contour lines are below or equal to the O+NO 
dissociation limit. The crossing seam is represented by the 
thick red line. The PES are represented as a function of the 
bending angle and the NO distance Ri, with the second NO 
distance R2 fixed to 1.19 A. The snapshots are obtained by 
representing only two coordinates (bending angle and Ri) and 
ignoring the third one (i?2). 



what simpler. The trajectories start on the excited adi- 
abatic state. Around t=10 fs, most of the trajectories 
jump on the ground adiabatic state and then remain in 
that state, in which a large amplitude bending motion 
occurs subsequently. Thus, the same dynamics is seen as 
a variation of the population in the diabatic representa- 



tion or as a variation of the average bending angle in the 
adiabatic representation. 

The time dependencies of the diabatic populations 
(Fig. [3fa) ) as well as the snapshots of the trajectories 
(Fig l4J) will be used in section IV to analyze the results 
of the experiment presented in the following section. 

III. HIGH HARMONIC TRANSIENT 
GRATING SPECTROSCOPY 

Transient grating spectroscopy is widely used in non- 
linear optical spectroscopy to perform background-free 
measurements. In this configuration, the molecular 
medium is excited by two synchronized non-colinear 
beams which optically interfere, forming a grating of 
molecular excitation. This technique can be directly 
transposed to the extremely non-linear regime of high- 
order harmonic spectroscopy^. The probe beam gen- 
erates high harmonics in this structured medium. The 
harmonics are diffracted by the grating of excitation, 
leading to the appearance of one or several diffraction 
peaks around the main harmonic beam. In the present 
experiment, the fringes of the excitation grating being 
horizontal the diffraction peaks appear above and below 
the main harmonic beam. This diffraction pattern can 
be used to perform high contrast pump-probe measure- 
ments, since there is no signal in the direction of the 
diffracted peaks in the absence of molecular excitation. 
In addition, since the diffraction pattern is formed by 
the interference of emission from excited and unexcited 
sources, it encodes the phase difference between ha rmon ic 
emission from excited and ground state molecule d 15 * 45 *. 

High harmonic transient grating spectroscopy is partic- 
ularly valuable in a context where the pump can induce 
dynamics both through one-photon and multi-photon ab- 
sorption. The spatial modulation of the pump beam in- 
tensity is sinusoidal. For a pure one-photon process, the 
excitation grating is also sinusoidal and only first order 
diffraction can be observed in the far field (Fig. |5ja)). 
If the excitation results from a two-photon process, then 
the grating is determined by the square of the sinusoidal 
modulation of the pump beam, which contains the sec- 
ond harmonic of the exciting beam spatial modulation 
frequency. Consequently, second order diffraction peaks 
appear (Fig. (5jb)). Note that second (and higher) order 
diffraction peaks can also appear in the case of a satu- 
rated one-photon transition: the saturation of the molec- 
ular excitation induces an anharmonicity of the grating 
which is reflected in the far field diffraction pattern (Fig. 
§b)). 

A. Experimental setup 

The experiment was performed using the Aurore Ti:Sa 
source at CELIA, delivering 7 mJ, 27 fs, 800 nm pulses 
at 1 kHz. The laser beam is split into a pump (20% 
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in energy) and probe (80% in energy) beams. The pump 
beam is frequency doubled using a 200 /im SH-BBO crys- 
tal, which provides wavelength tunability around 400 nm 
without reducing the temporal resolution (~ 35 fs cor- 
responding to a 7 nm bandwidth). The pump is fur- 
ther split in two parts by a 50-50 beamsplitter. The 800 
nm light is filtered out from the pump beams using two 
dichroic dielectric mirrors (in each arm) with high reflec- 
tivity around 400 nm and high transmission around 800 
nm. An adjustable aperture on the pump beams is used 
to reduce their energy and to ensure that the excitation 
of the molecular medium will be uniform over the area 
probed by high-harmonic generation. The pump energy 
used in the experiment typically varied between 10 and 
35 £tJ/p in each arm. The two pump beams are parallel 
and vertically shifted with respect to each other by 16 
mm. The 800 nm probe beam is centered between the 
two pumps. The three beams are focused onto the gas 
target using 750 mm radius spherical silver mirror placed 
under vacuum in order to minimize self phase modula- 
tion in the entrance window of the vacuum chamber. The 
pump beams produce an interference pattern at focus 
with a fringe spacing of 9.4 fim. The probe pulse waist is 
~35 /xm. The linear polarization of the probe pulse can 
be continuously rotated relative to the S-polarized pump 
pulses. The typical probe intensity determined from the 
harmonic cut-off (see Fig[5| is ~1.2xl0 14 W/cm 2 . The 
high-harmonic emission is produced a few 100 /im down- 
stream from a 100 /im nozzle gas jet (General Valve) 
operated at 100 Hz with a backing pressure of pure NO2 
at 780 mbar. To avoid the dimer formation, the nozzle 
temperature is fixed at 80° C where only ~ 10 % species 
are N2O4. The latter has an extremely low excitation 
probability at 400 nm (< 4%)P3, so that it does not con- 
tribute to the detected signal. The harmonic signal is 
dispersed horizontally by a grazing incidence spherical 
grating with variable groove spacing and imaged by a 
detector consisting of microchannel plates, a phosphor 
screen and a charge-coupled device camera. 



B. Experimental results 

As a first step we studied the picosecond dynamics 
of the diffracted and undiffracted signal. Similarly to 
refPHl, we observed picosecond transients reflecting the 
one-photon dynamics of NO2 with a 1.8 ± 0.3 ps disso- 
ciation time. This indicates that the dominant process 
detected in our measurement is dictated by one-photon 
excitation. 

In order to study the electronic relaxation of photoex- 
cited NO2 molecules shown in Fig. |3]and|4j we have per- 
formed pump-probe measurements in the first few hun- 
dred femtoseconds following excitation. The pump pulses 
were 20 /iJ each and centered at 400 nm. First we deter- 
mine angle between pump and probe polarizations that 
optimizes the contrast. Figure [6] shows the first order 
diffraction efficiency rj for harmonic 15 as a function of 
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FIG. 5. Principle of high harmonic transient grating spec- 
troscopy. In the case of one-photon excitation (a), the si- 
nusoidal modulation of the pump intensity in the near field 
produces a sinusoidal modulation of the harmonic emission, 
which results in the appearance of first order diffracted light 
up and down the harmonics in the far field. The far field pro- 
file was experimentally recorded by using two pump pulses 
containing 20 /iJ/p, at a pump-probe delay of 100 fs and with 
a probe pulse polarized orthogonally to the pump pulses. At 
higher pump intensity (35 p3/p, (b)), second order diffrac- 
tion peaks appear in the far field. They can result from two- 
photon transitions or from the saturation of the one-photon 
absorption, the two effects leading to anharmonicities of the 
near-field harmonic spatial modulation. 



delay and pump-probe polarization angle. The diffrac- 
tion efficiency is defined as the fraction of the harmonic 
light contained within a single diffraction order : 



V = 



Jo + I-i + h 



(4) 



Monitoring this quantity rather than the raw diffracted 
signal enables lowering the effect of laser fluctuations 
and provides smoother signals. The diffraction efficiency 
shows a maximum around 120 fs. The contrast of this 
feature is maximized when the pump and probe polariza- 
tions are perpendicular to each other. All of the following 
discussion will therefore be with reference to this config- 
uration. 

Figure [7] shows the evolution of the total signal (blue) 
and first order diffraction efficiency (red) as a function 
of the pump-probe delay (10 fs steps), for harmonics 15 
to 21. This result was obtained by adding up 5 consec- 
utive scans, resulting in an accumulation over 6250 laser 
shots per delay. The zero pump-probe delay is deter- 
mined from the even harmonics signal. These harmonics 
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FIG. 6. First order diffraction efficiency for H15 as a function 
of the angle between pump and probe polarizations. 
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FIG. 7. Time dependency of the total harmonic signal (blue) 
and first order diffraction efficiency (red, mutiplied by 20) for 
harmonics 15 (a), 17 (b), 19 (c) and 21 (d). The pump pulses 
are 20 /j.3 each, centered at 400 nm and polarized orthogonally 
to the probe. The cross correlation signal (panel (a): dots for 
the experimental points and green line for the gaussian fit) 
from harmonic 16 allows an accurate determination of the 
zero delay. 



result from wave mixing between the pump (400 nm) 
and probe (800 nm) pulses, and thus provide a cross- 
correlation signal. The influence of the dynamics initi- 
ated by the pump pulse on even harmonics is taken into 
account to determine the zero delay. We assume that 
the even and odd harmonics are similarly affected by the 
change of the molecular geometry. The cross correlation 
signal is thus obtained by normalizing the signal of even 
harmonics by that of adjacent odds. The resulting signal 
has a symmetric shape, which can be fitted by a gaus- 
sian function to extract the zero delay. Note that this 
normalization shifts the zero delay by +10 fs. 

All harmonics experience a fast intensity decay within 
the first few tens of femtoseconds of interaction with the 



pump. After t = 20 fs, the signal oscillates, presenting a 
maximum at 65 fs and a minimum at 130 fs. After that, it 
shows a slow increase, characteristic of the ongoing disso- 
ciation of a fraction of the excited molecules, since at 400 
nm, with a 7 nm laser pump bandwidth, 25% of excited 
molecules lie above the dissociation limit. The diffrac- 
tion efficiency shows an inverted behavior, with maxima 
at 15 fs and 120 fs, and a minimum at 58 fs. The con- 
trast of diffraction efficiency is better than that of the 
total signal, in particular for highest harmonics. In par- 
allel to the harmonic signal measurements, we measured 
the ionization yield from the generation jet (not shown 
here). The ion signal increases by a factor 1.3 between 
t = —60 fs and 60 fs and is stable afterwards. 

Overall, these results are similar to the previously re- 
ported observation^, nonetheless they present a few dif- 
ferences. First, the diffraction efficiency measured here is 
6 times lower than in the previous experiment. Second, 
dynamics appear as well in the undiffracted harmonic sig- 
nal and not only in the diffracted light. Third, only two 
oscillations of the diffracted light are observed: a first 
sharp peak, and a second one with a steep leading edge 
and slow falling edge. Last but not least, the exact tim- 
ing of the oscillations is different. In order to investigate 
the possible origins of these differences, we conducted a 
systematic study of the transient. 

When the pump energy is increased to 30 /iJ, sec- 
ond order diffraction peaks appear around the harmon- 
ics. Figure [8] compares the evolution of the first and sec- 
ond order diffraction efficiencies as a function of pump- 
probe delay. The dynamics observed for the second or- 
der diffraction is clearly the same as in the first order. 
This indicates that the second order diffraction peak orig- 
inates from the saturation of the excitation rather than 
from a two-photon absorption process. Indeed, with a 
2 B2 2 Ai absorption cross-sectiorP^l ~ 6xl0~ 19 cm 2 
and a pump waist size of ~ 60 /im, saturation is expected 
as soon as the pump energy reaches ^25 per arm. For 
the first order, the results are very similar to those ob- 
tained at lower pump intensity in Fig. [7] While the 
relative heights of the peaks at 20 fs and 130 fs slightly 
differ between the two measurements, their positions are 
identical. This shows that increasing the pump intensity 
does not affect the detected dynamics. We performed 
measurements varying the pump energy between 10 /xJ 
and 35 /iJ and did not observe significant changes in the 
temporal behavior of the first order diffracted signal. The 
first peak in the second order diffraction is slightly shifted 
to longer delays (Fig. |8jd)), which could reflect the dy- 
namics of the saturation process. 

As the rotational temperature of the molecules may 
play a role in the measurement, in particular through 
the dyna mics of the anisotropy induced by the pump 
excitatio n 1 48 ! 49 !, or the rovibrational and rovibronic cou- 
plings which are especially significant at long distances, 
we have compared the transients measured with two dif- 
ferent gas sources. The first source is that used for all 
results presented so far: a 100 Hz General Valve pulsed 
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FIG. 8. Time dependency of the first order (light red) and 
second order (dark red, mutiplied by 10) diffraction efficiency 
for harmonics 15 (a), 17 (b), 19 (c) and 21 (d). The pump 
pulses are 35 /iJ each, are centered at 400 nm and are polar- 
ized orthogonally to the probe. 
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FIG. 9. First order diffraction efficiency in various excitation 
conditions, (a) Signal obtained using a warm (red) and cold 
(yellow) molecular beam, (b) Influence of the pump wave- 
length: 400 nm (red), 395 nm (blue) and 392 nm (purple). 
The signals are vertically shifted for sake of clarity. 



source with a 100 /xm nozzle diameter, backed by ~ 780 
mbar of pure NO2, and heated at 80°C. In this configu- 
ration the laser is focused at a few hundred microns away 
from the nozzle. The second source is a 1 kHz Even-Lavie 
valve with a 250 /im conical nozzle, backed by a mixture 
of 16 bars of He and ~ 800 mbar of NO2, and heated 
at 120°C. In that configuration the laser is focused at 
~ 2 mm away from the nozzle. We expect the rotational 
temperature to be significantly lower in the latter case. 
Figure |9f a) compares the evolution of the diffraction ef- 
ficiency for harmonic 17 using these two sources. While 
the ratio between the first and second maximum is differ- 
ent, which probably reflects slightly different laser pump 
conditions, the results are remarkably similar. 

Last, we studied the evolution of the transient as a 
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FIG. 10. Evolution of the calculated bending trajectory 
packet for the diabatic surfaces 1 (a) and 2 (b), in logarithmic 
scale, (c) Comparison between calculated (dark) and mea- 
sured (light) high harmonic total signal (blue) and first order 
diffraction efficiency (red, multiplied by 20) for harmonic 15. 



function of the pump wavelength, by varying it between 
392 and 400 nm (Fig. ^b)). The results are robust 
against such wavelength variations. 

In summary, none of the external parameters like pump 
intensity, wavelength and rotational temperature, can ex- 
plain the differences observed between the data shown 
here and the ones published in refP^l. One possibility 
to explain the absence of oscillatory structure in the un- 
diffracted HHG signal and more persistent oscillations 
seen in the earlier work could be the noise level, the statis- 
tics being improved in the present work due to the higher 
repetition rate of the laser system. 



IV. INTERPRETATION 

The experimental results show that the evolution of 
the diffraction efficiency presents robust characteristics: 
a maximum around 20 fs, a minimum around 60 fs, and 
a broader maximum around 120 fs. This behavior can be 
compared to that of the molecular dynamics described 
in Fig. [3j The minimum diffraction efficiency at 60 
fs corresponds to a minimum in the population of the 
excited diabatic state, as well as to a maximum aver- 
age bond angle of the molecule. In order to distinguish 
which of the two plays a dominant role, we performed 
simple calculations of the harmonic signal and diffrac- 
tion efficiency. Since the tunnel ionization probabilities 
and recombination dipole moments have been calculated 
only as a function of the bending angle^, our model is 
one-dimensional: we only consider the influence of the 
motion along the bending coordinate by integrating the 
trajectory swarms over the stretching coordinates. The 
resulting evolutions are shown in Fig |10[ a-b). 

In a three-step picture of HHG, the harmonic signal 
is determined by the strong field ionization rate of the 
molecule, the propagation of the electron in the contin- 
uum and the recombination dipole moment of this elec- 
tron to the ion core. The ionization rates and complex 
recombination dipole moments were calculated in ReP^ 
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FIG. 11. Harmonic phase shift AI p t in units of 71", as a func- 
tion of bending angle, for emission from PES 1 (I 2 A') (blue) 
and PES 2 (2 2 A')(red). The shaded areas correspond to the 
main part explored by the trajectories. 



for the fundamental and first excited diabatic states, for 
three bond angles of the NO2 molecules (85°, 102° and 
134°), taking into account the different possible ioniza- 
tion channels. We use these results in our model and 
interpolate them to intermediate bond angles. The prop- 
agation of the electron in the continuum induces an addi- 
tional phase shift in the harmonic emission from excited 
molecules with respect to ground state molecules, which 
corresponds to the phase <p accumulated by the molecule 
between ionization and recombination : ip = A/p7^1£l] 
where AI p is the difference between the ionization po- 
tential of the considered state and the unexcited ground 
state, and r is the electron travel time, typically 1.2 fs 
for harmonic 17 for the experimental conditions. As the 
bond angle varies, the ionization potential changes be- 
cause of the dramatic difference between the neutral and 
cationic potential surfaces. This can induce important 
variations of the harmonic phase. We calculate this phase 
by using the ionic potential curves from ReP^. As shown 
in Fig. |11[ the phases vary linearly as a function of an- 
gle within the narrow range explored by the wavepacket, 
namely 120-142° for the ground electronic state and 90- 
110° for the excited state (shaded areas in Fig. 11). 

For each pump-probe delay i, the harmonic signal is 
calculated by coherently summing the contributions from 
the different bending angles 8, for each molecular state. 
The undiffracted and first order diffracted signals are 
given bjJ2H 



/ (i)cx J (l-r)d g + rW x (6,t)d x (d) 
+rW 2 (8,t)d 2 (0)d6 
Ji(*) oc i / ' rW x (9,t)(d x (9) - d g 

2 

rW 2 (9,t)(d 2 (9) - d g )d9 



(5) 



(6) 



where r is the fraction of excited molecules, d g the har- 
monic dipole moment for the non excited molecules, dj (9) 
the complex dipole moment for emission from surface 
j for the excited molecules, and Wj(0,i) is the time- 




100 200 
Delay (fs) 



FIG. 12. Total harmonic signal (blue) and first order diffrac- 
tion efficiency (red, multiplied by 20) obtained by imposing 
constant equal populations for the two diabatic states. The 
black dotted curves are the results of the full calculation. 



dependent bending trajectory packet on the j surface 
as depicted in Fig. 
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(a-b). Note that since the 
high harmonic signal is a coherent macroscopic emission, 
its amplitude is quadratically dependent on the popu- 
lation. The total signal and diffraction efficiency are 
hot = Jo + 2 x I x and rj x = J\jJ tat . Figure [Tf^c) shows 
the results obtained by assuming that r = 20% of the 
molecules are excited, together with the experimental 
results. The results of the calculations have been con- 
voluted by a gaussian function of 15 fs FWHM to take 
into account the duration of the harmonic generation pro- 
cess. The total signal is normalized to 1 at negative de- 
lays. The agreement is remarkably good, especially for 
the contrast, considering the simplicity of the model used 
for HHG calculations and the fact that we neglect the 
stretching motion. Both the fast dynamics and the av- 
erage level of signals at long delays are well reproduced. 
Besides the slight temporal shift in the position of the 
first peak in the diffraction efficiency, the main difference 
between experiment and theory lies in the width of the 
second oscillation: the experimental results show a single 
broad secondary peak around 120 fs, associated to a dip 
in the total signal, while the calculated signal shows a 
peak at 100 fs followed by weaker replicas at 170 fs and 
240 fs. This is probably an effect of the spreading of the 
packet along the stretching coordinates, which smoothes 
out the evolution at longer delays. 

In order to ascertain the origin of the transients ob- 
served on diffracted signals, we cancel out the effect 
of the population transfers by normalizing the pack- 
ets at each delay: P x (t) = fW x (9,t)d9 = 1/2 and 
P 2 (t) = fW 2 (0,t)d6 = 1/2. This enables us to isolate 
the influence of the bending motion. The results are re- 
markably similar to those obtained with the full calcu- 



lation (Fig. 12 1. The influence of the packet dynamics 



can be further disentangled by calculating the contribu- 
tion of each packet alone (Fig. 13 (a)), setting P x (t) = 1 
and P 2 {t) = (blue) or P x (t) = and P 2 {t) = 1 (red). 
This reveals that the main origin of the observed exper- 
imental modulations is the bending motion of the tra- 
jectory packet on the excited diabatic surface 2. This is 
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surprizing, given the fact that the extension of the bend- 
ing motion is more important on the diabatic surface 1 
than on 2. As the packets move, the harmonic emission 
is mainly modulated by the phase accumulated by the 
molecules during the electron travel in the continuum, 
AI p t (Fig. 11). This phase, which appears in the com- 
plex dipole moments, typically varies between — 0.4-7T and 
0.47T during the packet motion on surface 1, and between 
0.27T and 0.87T during the packet motion on surface 2. 
Figure 13 [b) shows the evolution of the diffraction effi- 



ciency as a function of AI p r. It is clear from this graph 
that the packet motion on surface 2 will induce a much 
stronger diffraction signal as well as a deeper modulation. 
This shows that similar variations of AI p t as a function 
of the reaction coordinates can have very different im- 
pact on the signal depending on the absolute value of the 
phase. For instance an additional recombination phase 
in the harmonic emission could strongly modify the situ- 
ation and make the dynamics on surface 1 dominate the 
signal. This issue is for instance illustrated in the dif- 
ferent contrast observed as a function of the harmonic 
order for the first and second order diffraction efficiency. 
The sensitivity of high harmonic spectroscopy to specific 
dynamics will thus depend on the considered system and 
must be carefully studied case by case. 

Since the signal is dominated by the vibrational motion 
on the excited diabatic surface, one may wonder what is 
the influence of the conical intersection on this signal. To 
estimate this, we calculated the molecular dynamics on 
surface 2 in the absence of electronic coupling to surface 
1. The results are shown in Fig. [T4|[b) and compared to 
the results including the coupling but keeping the pop- 
ulation constant P%{t) — 1 (FigpT^a)). The results are 
similar until 20 fs delay. This means that the first pas- 
sage through the conical intersection, around 10 fs delay, 
has only little effect on the structure of the trajectory 
packet. By contrast at 60 fs there is a strong transfer to 
the ground diabatic state (Fig [4]) which alters the packet, 
depleting the highest bending angles. This is followed by 
a transfer back to the excited state around 100 fs. Fig- 
c) compares the high harmonic diffraction efficicn- 
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cies obtained with (red) and without (blue) the coupling. 
The results are quite different. In particular they indi- 
cate that the second oscillation observed around 100 fs in 
the experimental signal is due to the transfer through the 
conical intersection from the ground to the first excited 
diabatic state. While the change in the overall popula- 
tion of the excited state has only little effect on the signal 
(Fig. [12] ), the packet structure and dynamical changes 
induced by transfers through the conical intersection are 
clearly reflected in the signal. 

The present results have to be compared to recent 
calculations based on wavepacket propagation^. In a 
first ID calculation, the bending wavepackets were mod- 
elled as Gaussian functions of constant width, with their 
centers following the angular position of the full 3D 
wavepackets center. The harmonic emission was calcu- 
lated using the same procedure as described here, i.e. 



taking into account the phase modulations induced by 
the vertical ionization potential changes, and the changes 
in the harmonic dipole moments due to the change in 
molecular geometry. This calculation resulted in strongly 
contrasted modulations of the diffracted signal as a func- 
tion of delay. In a second calculation, the wavepacket 
motion was neglected and only the influence of the pop- 
ulations was retained. The resulting modulation of the 
diffracted light was weaker than in the ID case. Last, a 
3D calculation taking into account the evolution of the 
molecular wavepackets along the three dimensions was 
performed, assuming constant harmonic dipole moments 
and phases except for from the I p t phase. The results 
hardly show any modulation in the harmonic signaP^. 
Since the model neglecting wavepacket motion and tak- 
ing into account only the population evolutions gave the 
best qualitative agreement with their experiments, Kraus 
et al. concluded that the experimental signal reflected 
the variation of the diabatic state populations, the phase 
modulations I p t associated to the vibrational motion be- 
ing quickly washed out by wavepacket spreading in the 
three dimensions. This conclusion is quite different from 
that drawn from our ID model, in which we find that the 
population transfer only has a weak effect on the detected 
signal. It is interesting to note that for the ID models, 
the proper treatment of the extension of the wavepacket 
along the bending angle (as done in equation 9 and 10 
in the present paper) without the assumption of a Gaus- 
sian shape as done irP^ leads to a better agreement with 
experiments. 

If only the populations or vertical ionization potential 
were probed by high harmonic spectroscopy, the mea- 
surement should reveal the same dynamics irrespective 
to the pump-probe angle. This is not the case: the dy- 
namics is much more contrasted for orthogonal pump and 
probes, as shown in Fig. [6j This is a consequence of the 
angular dependence of the harmonic signal with respect 
to the z- molecular axis: the harmonic dipole moments 
can indeed play a major role in the outcome of the mea- 
surement. The fact that the ID model used here repro- 
duces better the experiment than the 3D model previ- 
ously used irP^l assuming constant harmonic dipole mo- 
ments for all bending angles and bond lengths suggests 
that the harmonic dipole moments present significant 
variations along the stretching coordinates, such that the 
complete 3D averaging is closer to the complete ID cal- 
culation than to the partial 3D calculation. A complete 
3D calculation, including the evolution of the harmonic 
dipole moments along all molecular coordinates, is re- 
quired to confirm this statement. 



V. CONCLUSION 

High harmonic transient grating spectroscopy experi- 
ments in photoexcited NO2 molecules show modulations 
of both the total harmonic signal and first order diffrac- 
tion efficiency within the first 200 fs following the pump 
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FIG. 13. (a) Partial contributions to the diffraction efficiency, 
considering only the trajectories on surface 1 (Pi(t) = 1, 
P 2 (t) = 0, blue) or surface 2 (Pi(t) = 0, P 2 (t) = l,red). 
The black dashed line is the result of the full calculation, (b) 
Normalized diffraction efficiency as a function of the phase 
modulation AI p t in units of n. The shaded areas correspond 
to the main part explored by the trajectories. 
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FIG. 14. Temporal evolution of the calculated trajectory 
packet on the excited diabatic surface 2 with (a) and without 
(b) the coupling to the ground state, (c) First order diffrac- 
tion efficiency induced by the trajectories on surface 2, with 
(red) and without (blue) the coupling to the ground state. 
The black dashed line is the result of the full calculation. 



pulse. These can be remarkably well reproduced by com- 
bining Trajectory Surface Hopping calculations of the re- 
action dynamics around the conical intersection and a 
simple model of high harmonic generation. The anal- 
ysis reveals that the dominant information encoded in 
the transient grating signal is the molecular vibration on 
the first excited diabatic state. This sensitivity results 
from the modulation of the vertical ionization potential 
as the molecule vibrates, which induces a phase modula- 
tion of the harmonic emission. The vibrational motion in 
the ground diabatic state also induces a significant phase 
modulation. However this modulation is centered on zero 
(the zero phase shift corresponding to the equilibrium ge- 
ometry of the molecule in the ground state), such that it 
has only little effect on the measured signal. 

These results show the importance of the harmonic 
phase in the nature of the information encoded in high 
harmonic transient grating spectroscopy experiments. 
Indeed, an additional phase shift in the harmonic emis- 
sion could invert the situation and make the signal most 



sensitive to the vibration in the ground diabatic state. 
Since the phase of the recombination dipole moments in 
HHG can s how si gnificant structures as a function of har- 
monic ordei^^, it would be interesting to extend this 
study to longer driving wavelengths which produce much 
broader harmonic spectra. The transient grating signal 
could then encode different dynamics (ground diabatic 
state vibrations, first excited diabatic state vibrations, 
or population transfer through the conical intersection) 
depending on the harmonic order, giving a full picture of 
the ongoing dynamics. 
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